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Abstract 

Nucleon form factors are calculated in a relativistic diquark-quark picture based 
on the Nambu-Jona-Lasinio model. The nucleon wave function is obtained in a 
static approximation to the quark exchange interaction between the valence quark 
and the diquark. We evaluate the valence quark and 0 + -diquark contribution to 
the nucleon electromagnetic and weak currents. We find reasonable electric charge 
radii, magnetic moments as in the additive diquark model, and gA ~ 1. We discuss 
the dependence on the model parameters. 
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A quantitative understanding of baryon form factors at low energies in terms of quark 
degrees of freedom remains a challenging problem. Attempts to base a low-energy de¬ 
scription of baryons on QCD usually rely on the limit Nq —> oo, in which a soliton picture 
emerges. While properly describing the effects of the non-perturbative meson cloud, this 
approach suffers from difficulties related to the breaking of translational invariance. For 
many purposes, such as e.g. the calculation of timelike nucleon form factors, it is im¬ 
portant to maintain proper relativistic kinematics. A description of the nucleon as a 
relativistic bound state can be realized in quark models with an effective 2-body inter¬ 
action 1B.B- For Nq = 3, this naturally leads to the concept of diquarks as effective 
constituents |4|]. The simplest such model is a Nambu-Jona-Lasinio (NJL-) model with a 
local (current-current) interaction in the 3-color channel. A corresponding color-singlet 
interaction describes the spontaneous breaking of chiral symmetry and provides a reason¬ 
able account of meson properties Hi- The NJL model can be rewritten as an effective 
theory of mesons and baryons ||, in which bound states of composite diquarks and quarks 
occur due to quark exchange. The spectrum of | + -baryons has been calculated using a 
static approximation for the quark exchange interaction |7|]. The Faddeev equation for 
the nucleon has also been solved exactly P , including also the l + -diquark channel |9], [10[ . 


Here, we apply this description to the study of the nucleon form factors at low momen¬ 
tum. We consider the simplest case of a bound state of a 0 + -diquark and a valence quark 
and take the baryon wave functions in the static approximation to the quark exchange. 
Given the complexity of any relativistic description of baryons, it is useful to first explore 
this simple approximation. The present study extends a previous analysis of the form 
factors in an additive diquark-quark scheme [[□]]. Contrary to that approach, we now 
derive the relativistic nucleon wave function within the model quark dynamics in a crude 
but consistent approximation. Our intention is to estimate the role of diquark and valence 
quark currents in the nucleon electromagnetic and weak form factors. Moreover, we wish 
to point out some general features of the bound state description of baryons, which are 
likely to persist in more elaborate calculations with exact Faddeev wave functions. 

The basis of our approach is the effective baryon action, which is obtained by rewriting 
the partition function of the NJL-model in terms of composite diquark and baryon fields 
0, [T. It is defined in terms of the baryon Greens function, 


G b = [Go 1 -H}~ + a G b . 


( 1 ) 


Here, 

G 0 (k,k';P) = V(\P - k)G{\P + k')5(k - k') (2) 

is the product of the scalar diquark propagator, T >, which is the solution of the diquark 
Bethe-Salpeter equation of the NJL model 0, and the quark propagator, G[k ) — — 
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M) 1 . We consider the isospin limit, M u = Md = M. P and k are the total and relative 
four-momentum of the system. Furthermore, 

H(k, k') = YG(—k — k')Y (3) 


is the quark exchange interaction. For a scalar diquark, F = iqsCe, where C is the charge 
conjugation matrix and e the generator of the 3-representation of the color group. The 
contact term in eq.(l|), A G B , will not be important in the following. Baryon masses and 
bound state wave functions are obtained as color-singlet solutions of the Faddeev-type 
equation 

d A k' 

G^ik.k'-P^P) = 0, P 2 = m 2 B . (4) 

In the static approximation for the quark exchange 10], 

H(k,k') - -iA r , (5) 

the baryon wave function in eq.(f|) can be written as 

ip(k\P) = V(\P -k)G{\P + k)<t>(P). (6) 


The baryon amplitude, 4>(P), is then obtained as the solution of a Dirac-type equation, 

( A(P 2 )If + B(P 2 ) - M) <f>(P) = 0, (7) 


where the effective Dirac operator is defined as 

r d 4 k 

A(P 2 )If + B(P 2 ) = / V(\P-k)G(\P + k). (8) 

In this approximation the baryon amplitude does not depend on the relative momentum 
of the diquark-quark pair. The integral in eq.(0) is defined with a momentum cutoff, A#, 
which effectively controls the width of the baryon wave function. (A# is not identical to 
the NJL cutoff, cf. below.) We normalize the on-shell baryon amplitude by the condition 


}(p) 


pv 


d 

dP„ 


(A(P 2 )I? +B{P 2 )) 


4>{P) = m B - 


(9) 


Eq. (0) fixes the correct relativistic dispersion law for the c.o.nr. motion of the bound state 
and is identical to the usual normalization condition for Bethe-Salpeter amplitudes. 

Baryon form factors are calculated by coupling to the free baryon action an electro¬ 
magnetic or weak gauge held, by minimal substitution at quark level. One then 
obtains the current of the baryon held as the corresponding functional derivative of the 
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inverse baryon propagator, eq.(fl). The matrix element of the baryon electromagnetic 
current for a transition P ; —> Pf, where P 2 = P f 2 = m 2 B , is in general given by 

c d A k r d A k' ~ 

= ( n(Q 2 ) - p B (Q 2 )) 2 m B |hit + Fg(Q 2 )Y. (11) 

Here, Q = P{ — P ■ and F B ,F B are the baryon electric and magnetic form factor. The 
baryon current consists of terms describing the coupling of the external field to the diquark 
propagator, the quark propagator and the quark exchange. (The contact term in eq.([T|) 
does not contribute to the baryon current on the mass shell.) In the static approximation, 
eq.([5]), it is consistent to neglect the quark exchange current. Note that the quark exchange 
also does not contribute to the normalization of the baryon amplitude, cf. eq.(^). Thus, 
in the static approximation, eq. (0) becomes 


r B (p<, p ,) = m) (r wn, p i) + r B,„» rk(n, a) + c(a b ) 7 «) ^ (12) 

r d 4 k 

r^ diq (Pf, Pi) = J AbWT ^ - - fc, |Pi - k)v(ln - k)G(\Pt + iPi + k) 

d A k 


I * rj^ H 

p 0 = J AB -^y v ^ P{ + i Pi - fc ) G ^ Pf + 


(13) 

(14) 


We note that the on-shell baryon current defined by eqs.(|l|, |lj) is transverse. Here, 


F^k, k’) = ( Q u + Q d )Z~ A {{k + kyF 0+ (q 2 , k 2 , k’ 2 ) + q^G 0+ (q 2 , k 2 , k' 2 )} , (15) 


with q = k—k 1 , is the off-shell e.m. vertex function of the extended 0 + -'ud-diquark HTTjp^]. 
The valence quark is taken as pointlike, T^ = iQ utd r y tl , corresponding, respectively, to a 
proton or a neutron. (Q U}d = §, — | are the quark charges.) In the approach based on 
||, a valence quark form factor would arise if dynamical meson exchange were taken into 
account^. Furthermore, we have included in eq.(|T^) a counter term, C(A B )'y IJ ', which 
we choose to fix charge conservation at Q 2 = 0. This is necessary when working with a 
momentum-space cutoff. However, in the numerical calculations we find that the violation 
of charge conservation is small, he., C(A B ) -C 1. Note also that the counter term does 
not affect the charge radii. 

The composite diquark propagator and the corresponding e.m. vertex function, eq.([15|) 
have been determined within the gauge invariant proper-time regularized NJL model 


: A consistent treatment of meson exchange effects would require the inclusion of meson exchange in 
the quark-quark interaction and the diquark-quark binding 


4 











Tl| , p~2| . To facilitate the evaluation of the loop integrals in eqs.(|13|, [L4|), we approximate 
the full cliquark propagator by its pole form, 


Vik) 


Z 0 + 


k 2 — rn,Q + 


(16) 


where the pole and the residue are determined from the exact solution. This is justified 
if the cutoff defining the baryon wave function, A#, is not too large, which is the case 
in the calculations below. Consistently with the pole form of the propagator, eq. (H§, 
we make an on-shell approximation^ for the composite diquark vertex function, he., we 
replace in eq. © the diquark form factor F 0 +(q 2 , k 2 , k' 2 ) by its value at k 2 = k' 2 = 

(The longitudinal part of eq.(|T5D does not contribute to the on-shell baryon form factor.) 
For regularization of the diquark-quark loop integrals in eqs. ©0) we employ a sharp 
euclidean cutoff, A B , after introducing Feynman parameters, cf. ||. (We have also done 
calculations with a proper-time-type regularization for the diquark-quark loop and found 
only very small differences to cutoff regularization.) The contributions of eqs. ([Td|. [L4|) to 
the nucleon electric and magnetic form factor can then be identified after making use of 
the Dirac equation, eq.(|7|), to simplify the numerators of the Feynman integrals. 

Before discussing results for the baryon charge radii and magnetic moments, it is 
instructive to consider the dependence of these quantities on the baryon c.o.m. momentum. 
Fig.l shows, for a fixed quark and diquark mass, the proton charge radius from eq.([l2|) 
as a function of P 2 = P 2 = P 2 ^ m 2 B . Here, the amplitudes are normalized according to 
eq.(||), but at P 2 ^ m 2 B . (One may imagine the off-shell values of P 2 to be generated by 
varying the strength of the static quark exchange term, M, in eq.(^i).) The charge radius 
grows and eventually diverges, as P 2 approaches the diquark-quark continuum threshold. 
The presence of a singularity is, of course, due to the lack of confinement of this model. 
The strong variation of ( r 2 ) p with P 2 for bound baryons, however, should not be regarded 
as unphysical. We know that even in confining theories there are “would-be”-thresholds 
related to losely bound states [|T3|. In this spirit, fig.l simply reflects the fact that, in 
the context of the effective quark dynamics described by the NJL model, the baryon has 
to be considered as a losely bound state, contrary to e.g. the pion. This is in agreement 


with the fact that the diquark masses obtained in the NJL model [Q as well as in the 
instanton vacuum [Hj are typically of the order of 0.5... 0.6 GeV, so that the diquark- 
quark threshold is close to the physical nucleon mass. Furthermore, the exact Faddeev 
calculation leads to nucleon binding energies of e B = M + m 0 + —m B ~ 50 ... 100 MeV -C 
m B Q. When fixing the parameter A# of the static approximation, we must keep in mind 
that it is the binding energy — not the absolute mass — which is smoothly related to 

2 Note that the off-shell dependence of the exact V(k) and k') cancels to leading order in ( k 2 —m 2 d ) 
and (k 2 — m 2 ), when considering the product k')V{k') Q. 
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the characteristics of the baryon wave function. We thus take values of A B , for which the 
Dirac equation, eq.(^) gives binding energies comparable to those of the exact solutions 
(see fig. 1). These lead to realistic charge radii. We see no physical grounds for identifying 
A b with the NJL cutoff employed in the diquark and meson sector 0. the more since 
values of A^ ~ 0.6...0.7GeV would considerably overestimate the binding energy, cf. 
hg.l and 0. For a given value of A B , we then choose the NJL coupling constant in 
the diquark channel such as to reproduce the physical nucleon mass in eq.([7|) [[?[]. We 
emphasize that, in this way, the baryon form factors are not sensitive to the precise value 
of the baryon mass, if Ab (and thus, approximately, e B ) is kept fixed. 

With the parameters of the model determined in this way, we may now discuss the 
results for the on-shell nucleon charge radii and magnetic moments. As expected, the 
proton charge radius depends inversely on the cutoff used in the diquark-quark wave func¬ 
tion. The experimental radius, ( r 2 ) p = 0.74 fm 2 , can be reproduced for A B ~ 0.32 GeV. 


we 


Since this value for A B would probably be larger, if the l + -diquarks were included []TT 
give in table 1 also results for larger values of A B - In addition, we have performed calcula¬ 
tions with a pointlike diquark form factor, F 0 + = 1, cf. eq. Comparison of the values 
obtained with extended and pointlike scalar diquarks shows that the intrinsic charge ra¬ 
dius of the scalar diquark contributes about 10 ... 20% of the proton charge radius. The 
neutron charge radius is reproduced well. This quantity originally provided evidence for 
scalar diquark correlations [|T6| . The intrinsic diquark radius somewhat increases the neu¬ 
tron charge radius, improving the agreement with the experimental value. The magnetic 
moments are not too different from their values in an additive diquark-quark picture 
with only 0 + -diquarks, /i p = |(e/2 M),/i n = — A.(e/2 M). Here, large contributions are 
to be expected from the l + -diquark channel. Contrary to the charge radii, the magnetic 
moments are rather independent of the size of the baryon wave function. 

We have also calculated the isovector axial coupling constant of the nucleon^, g a ■ In 
this case the quark vertex is given by . The scalar diquark has no axial 

coupling, so that the baryon axial form factor is determined entirely by the corresponding 

now.) 


valence quark axial current, eq. (|14l) . (We do not include a counter term in eq. 

As shown in table 1, for realistic parameters, g^<l. The value Qa — 1 would be obtained 
in an additive diquark-quark scheme with only 0 + -diquarks. Corrections to this value 
are seen to increase with A B , as the diquark-quark system becomes more strongly bound. 
Note that also in we expect sizable contributions from the l + -diquark channel, both 
through the axial coupling of the l + -diquark as well as through 0 + -l + transitions. 


3 We do not address here questions related to the induced pseudoscalar form factor and the Goldberger- 
Treiman relation, which require the inclusion of meson fluctuations, but stay within a valence quark 


picture. See 0 for a discussion of these issues. 


6 







In summary, our results show two types of behavior of nucleon observables. Charge 
radii depend strongly on the extent of the diquark-quark wave function and can reasonably 
be described in a picture with (even pointlike) 0 + -diquarks. On the contrary, magnetic 
moments and (Ja are rather insensitive to the strength of the diquark-quark binding and 
therefore amenable to an additive cliquark-quark description JTT]|. They require, however, 
the inclusion of the l + -diquarks. 

In our simple estimates, we have treated the cutoff, A B , as a parameter and related 
it to the proton charge radius. In principle, the width of the diquark-quark wave func¬ 
tion should be determined by the underlying effective quark dynamics. It remains to be 
seen whether the quark exchange mechanism included in eq.(f|) leads by itself to realis¬ 
tic baryon radii, or whether these quantities are dominated by other effects (dynamical 
meson exchange, confinement). Meanwhile, it is of practical interest that, with appropri¬ 
ate values of A#, the static approximation seems to offer the possibility of a reasonable 
approximation to the nucleon wave function, if l + -diquarks are included. It may e.g. be 
employed to calculate strong form factors of the nucleon. 

The authors are grateful to H. Reinhardt and A. Buck for helpful discussions. 
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A b = 0.323 GeV 

A b = 0.35 GeV 

A b = 0.45 GeV 

exp. 

92/91 

1.46 

1.42 

1.27 


m 0 + / GeV 

0.582 

0.593 

0.631 


71/2 

z o+ 

14.3 

14.0 

13.0 


e B /GeV 

0.044 

0.055 

0.093 


(r 2 ) p / fm 2 

0.74 (0.66) 

0.65 (0.57) 

0.48 (0.40) 

0.74 

(r 2 ) n /fm 2 

-0.11 (-0.19) 

-0.10 (-0.17) 

-0.07 (-0.15) 

- 0.12 

H P /(e/2m B ) 

1.67 

1.67 

1.75 

2.79 

li n /(e/2m B ) 

-0.78 

-0.82 

-0.83 

-1.91 

9A 

0.96 

0.95 

0.91 

1.26 


Table 1: The nucleon electric charge radii, ( r 2 ) P)U , magnetic moments, n PtH , and isovector 
axial coupling constant, gA, for various values of the cutoff in the diquark-quark momen¬ 
tum, Ab- Results are given for a constituent quark mass of M = 0.4 GeV. For each value 
of Kb, the coupling constant of the scalar diquark channel, g 2 , has been chosen to fit the 
physical nucleon mass, tub = m p = 0.938 GeV. The corresponding values ofm 0 +,Z 0 +, 
and the diquark form factor have been determined from the full diquark propagator of the 
proper-time regularized NJL model 0 - The ratio of the diquark coupling constant to 
that in the meson channel, g 2 /gi, is as defined in [ 0 , 0 - Also given is the baryon binding 
energy, e B = M + m 0 + — m p . With a value of K B = 0.323 GeV, the experimental proton 
charge radius can be reproduced. The charge radii in parentheses have been calculated with 
a pointlike diquark form factor, F 0 + = 1, cf. eq. ( 0 ). 
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Figure caption 


Fig. 1 : The dependence of the proton charge radius, ( r 2 ) p on the proton c.o.m. momen¬ 


tum, P 2 = Pf = P{ 7 ^ rrf B . The charge radius is calculated from eq.(|l2D for arbitrary P 2 , 


with amplitudes normalized according to eq.@, but at P 2 7 ^ m B . Shown are the radii ob¬ 
tained with Kb = 0.55 GeV (solid line), 0.45 GeV (dashed line) and 0.35 GeV (dotted line). 
The on-shell masses corresponding to the different values of A B are marked by squares 
on the respective curves. Here, for all A B , the constituent quark mass is M = 0.4 GeV, 
and for the diquark mass we have chosen a fixed value of m<i = 0.6 GeV, so that the 
diquark-quark continuum starts at 1.0 GeV. 
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